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ABSTRACT 

Waves  in  the  surf  zone  are  a  highly  nonlinear  process 
which  is  evident  by  the  appearance  of  secondary  waves.   The 
secondary  waves  appear  as  strong  peaks  in  the  period  PDFs 
corresponding  to  the  first  harmonic  of  the  peak  of  the  wave 
spectrum.   The  strong  first  harmonic  period  peak  is  also 
reflected  in  the  highly  correlated  height  and  velocity  PDFs. 

Ereaking  waves  may  be  viewed  as  a  gradation  between  two 
extremes.   The  period,  height  and  velocity  PDFs  for  spilling 
breakers  generally  were  unimodal,  whereas  for  plunging 
breakers  strong  bimodality  was  found. 

The  joint  probability  density  functions  for  periods  and 
heights  of  the  breaking  waves  show  high  correlation  (0.60- 
0.80)  which  says  that  greater  wave  periods  are  associated 
with  larger  breaker  heights.   The  joint  PDFs  of  period  and 
particle  velocity,  and  velocity  and  height,  suggest  that  the 
maximum  onshore  particle  velocities  are  correlated  with  both 
the  wave  periods  and  wave  heights. 

Due  to  the  high  probability  of  the  secondary  waves,  the 
mean  wave  period  for  breakers  is  a  poor  descriptor  of  the 
average  period  of  the  offshore  incident  waves. 
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I.       INTRODUCTION 

Waves    in    deep  water   can   generally  be    represented   to   a 
first   approximation   as    a   Gaussian   process   having   zero   skew- 
ness    and  kurtosis    equal    to    three.      The   random  sea   surface 
is   viewed   as    the    linear   superposition   of   statistically   inde- 
pendent  spectral    components   which   is    adequately    described 
by   the   energy   density   spectrum.      The   restriction   of   statis- 
tical  independence    implies    the   wave    components    act    as    free 
waves    and   there    is   no    coupling   or   interaction   between 
various    spectral    components. 

As    the   waves   begin    to    shoal,    the    crests    become   more 
peaked,    and   the    troughs    become   more    elongated;    this    results 
in   a  positive   value    of   the    third   statistical   moment,    i.e., 
the   probability    density    distribution   of   sea   surface    eleva- 
tions  becomes    skewed   to    the    right.      The    skewness    is    a   result 
of  phase    coupling   of   the   wave    components,    which   is    indica- 
tive   of   the    increasing   importance    of   the   non-linearities 
of    the    system.       This    initially  weakly   non-linear   system 
found   in   shallow  water   may   be    adequately    described  by    a 
Gram-Charlier  probability    distribution    (Longuet-Higgins , 
1963). 

Near   breaking   waves    become    asymmetrical    about    the    crest 
with    a    continually    steepening   of    the   wave    face.      Killer 
(1976)    has    shown    that    a   jet    is    then    formed    at    the   wave 
crest    at   breaking,    which    extends    forward   and    then    down 


until   it    closes   with    the    forward   face   of   the   breaking  wave 
form.      The   asymmetry    about    the  wave    crest   yields    a    change 
in   the    fourth   moment,    or  kurtosis,    which   is    also    indicative 
of   the    increasing  non-linearities.      Thornton   et    al.     (1976) 
have  noted   that   secondary  wave    formation   becomes    apparent 
on    the    trailing   face   of   the  wave    form,    due    to    transfer   of 
energy    from  the   primary    frequency    to    its    harmonics   which 
also   points    to    a   strongly  nonlinear   system. 

This    thesis    uses    data   collected   as    described  by   Galvin 
(1975)    and  Olsen    (1977)    plus    additional    data   collected  by 
the    author.      The   experiments    were    designed   to   measure    the 
sea   surface   elevation    and   the   onshore   particle   velocities 
by  means    of   capacitance   wave    gauges    and   electromagnetic 
flowmeters,    respectively.      The   experiments   were    conducted 
at    a  variety   of  beaches    in   order    to   include   wave    types    of 
spilling,    plunging,    and   collapsing    forms.      The    characteris- 
tics   of   the   breaking  wave    forms    are    a    function   of    the   off- 
shore wave    spectrum   and   the  nearshore   bottom   slope    and 
topography,    where    the   increasing   beach   slope    increases    the 
rapidity    of    the    breaking   of    the   waves. 

Empirical  distributions  were  obtained  for  the  periods, 
wave  heights,  and  maximum  onshore  particle  velocities,  for 
waves  which  had  been  defined  by  the  zero-up-cross  method. 
Each  of  these  distributions  was  compared  with  the  Rayleigh 
probability  density  function  (PDF).  Joint  PDFs  of  periods 
and  wave  heights,  maximum  velocities  and  wave  heights,  and 
periods    and   maximum  velocities    were    calculated. 


II.   HISTORICAL  BACKGROUND 

A.   QUALITATIVE  DESCRIPTION 

Thornton  et  al.  (19  76)  have  described  a  number  of  ob- 
served universal  similarities  for  the  various  types  of 
breakers  occurring  on  different  beaches.   Normally  there  is 
a  quick  drawdown  of  water  just  before  the  breaker  arrives, 
followed  by  a  steep,  nearly  vertical  leading  edge.   Miller 
(1976)  suggests  that  subsequent  to  the  development  of  the 
leading  vertical  edge,  a  jet  forms  at  the  crest.   This  jet 
progressively  extends  itself  forward  and  down  until  it  unites 
with  the  forward  slope  of  the  breaking  wave,  resulting  in 
the  formation  of  a  vortex  with  a  circulation  around  a 
cavity.   When  the  tip  of  the  "jet"  touches  the  forward  slope 
of  the  breaking  wave,  the  resulting  force  forms  a  second 
weaker  vortex  in  front  of  the  first  one.   In  this  manner,  a 
series  of  vortices  yields  a  series  of  large,  single  pulse 
eddies,  which  contribute  to  the  rapid  decay  of  the  wave 
during  breaking.   It  has  been  found  that  the  magnitudes  of 
breaker  vortices  are  generally  correlated  with  a  gradation 
from  plunging  to  spilling  breakers. 

Thornton  et  al.  (1976)  have  observed  a  sloping  profile 
toward  the  trailing  edge  of  the  wave.   Secondary  waves  are 
frequently  noted  on  the  trailing  edge  both  on  the  analog 
strip  charts,  Figure  1,  and  by  visual  observations  in  the 
field.   The  secondary  waves  have  been  observed  to  grow  during 
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Figure  1.   Strip  Chart  Record  of  Sea  Surface  Elevations 
and  Particle  Velocities. 


11 


the  shoaling  process,  especially  just  prior  to  breaking 
when  the  primary  waves  are  steepening  significantly.   The 
secondary  waves  are  a  result  of  the  transfer  of  energy  from 
the  primary  wave  frequency  to  harmonics,  and  help  maintain 
the  potential  energy  across  the  surf  zone. 

Gallagher  (19  71)  has  photographed  secondary  wave  forma- 
tion occurring  across  reefs.   He  suggests  that  the  secondary 
waves  are  brought  about  by  a  response  to  vertical  accelera- 
tion near  the  front  of  a  broken  or  breaking  swell.   When  the 
amplitude  and  surface  curvature  of  a  shallow  water  wave  ex- 
ceed certain  values,  the  non-linearities  of  the  vertical 
accelerations  bring  soliton  wave  forms  into  being. 

Figure  1  shows  that  the  maximum  particle  velocity  occurs 
near  the  time  of  maximum  surface  elevation.   It  can  also  be 
noted  in  Figure  1  that  there  are  secondary  velocity  maxima 
which  correlate  with  the  secondary  surface  elevation  maxima. 
The  secondary  velocity  peaks  also  are  a  result  of  the  non- 
linearity  associated  with  the  breaking  wave. 

B.   EIIERGY  DENSITY  SPECTRA  OF  BREAKING  WAVES 

In  order  to  better  understand  the  PDFs  as  determined  in 
this  study,  the  energy  density  spectra  measured  by  Galvin 
(1975)  and  Olsen  (1977)  were  reviewed.   The  energy  density 
spectrum  is  proportional  to  potential  wave  energy  within  a 

frequency  band.  The  energy  density  spectrum,  G  (f) ,  indicates 

2 
how  the  variance,  a   ,  is  distributed  over  frequencies,  f, 

where 
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ax2   =   /         Gx(£)df  (1) 


(Kinsman,    1965), 

The   energy    density   spectra,    as    estimated   from  the    sea 
surface   elevation   records,    were    used   to    find   the    frequency 
of  peak    density   and   the    occurrence    of  harmonics.      The   power 
spectra   that  were  measured  by   Galvin   showed   a  peak    density 
over   the    frequency    range    of   0.061    to    0.092    Hz    (16.4    to    10.9 
sec).      First,    second   and   third  harmonics    can   be    observed 
as    shown   in    the   spectra   of   surface    elevation   and  horizontal 
velocities    for   4  March    19  75.      Additionally,    subharmonics 
were   observed  at    0.011    and   0.922   Hz    in    the    velocity   spectra. 
Galvin    (1975)    has    suggested   that    the    observed   low   frequency 
sub-harmonic  peaks   may   be    due    to    edge   waves    in    the    surf   zone. 

Thornton   et    al.     (1976)    note    that    the   appearance    of   strong 
harmonics    could  be    real    due    to    secondary  waves    or   artificial 
due    to   Fourier   analysis    of   a  peaked  wave    form.       It    is   noted 
that    secondary   maxima   of   n    are    apparent   in    the    strip    chart 
record   shown    in   Figure    1.       It    is    felt    that    the    secondary 
spectral   peaks    are    real    since    secondary  waves    appear    on    the 
strip    charts. 

G.       PROBABILITY   DENSITY   FUNCTIONS    OF   PERIODS    AND    HEIGHTS 

Studies    of    the    PDFs    of  period,    height,    and  wave    lengths 
have   been    conducted   over    the    last    20    years.       The   majority 
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Figure  2.   Energy  Density  Spectra  of  n  and  v,  4  March  19  75. 
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of   these    studies   have   been  mainly    concerned  with    theoretical 
or   deep   water   aspects    of    the   problem. 

The    theoretical    distribution   of  wave    amplitudes    for   a 
narrow-banded  Gaussian   surface    elevation   in    deep   water,    was 
shown  by   Longuet-Higgins    (1952)    to   be    a  Rayleigh    distribu- 
tion.     Tayfun    (1977),    in    studying    the    transformation   of   deep 
water  waves    to    shallow  water  waves,    showed   that    the   Rayleigh 
distribution    for  wave   amplitude   was    generally    applicable    to 
all   bandwidths.       Using    the    assumption    that    the   wave    height 
is    twice    the   wave    amplitude,    the   wave   height    PDF    is    then    re- 
presented by: 


p(H)dH   =    ^  exp 


8a 


dH  (2) 


where      a      is    the   variance    of    the   surface    elevation.      Using 
pressure    records    in    the    Gulf  of  Mexico,    Longuet-Higgins 
(1975)    observed   that    the   Rayleigh    distribution    fit    the    ob- 
served  distribution    reasonably  well.       He    found    that    there 
is    a   slight   excess    of  waves   with   heights   near   the   middle 
of   the    range    and   a    deficit    at    the    two    extremes.       Since    much 
of   the   high    frequency   portions    of    the  wave    records    are    fil- 
tered out   by    the   pressure    transducer,    Longuet-Higgins    (1975) 
suggests    that    the    narrow-band   approximation    may   not    be    as 
applicable    for    the    unfiltered   records.       In   shallow  water 
with   much    steeper  waves,    the    distribution    can    again   be 
expected   to   be    less    applicable    due    to    the    non-linearities 
which   become   increasingly    important. 
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Since    the   statistics    of     H      and     T      for   a  broad-banded 
spectrum  were    still   somewhat   unclear,    Goda    (1970)    simulated 
wave   profiles    using   an   infinite    series   where    the    amplitudes 
were    specified  by   various    theoretical    spectra   and    the   phase 
was    random  with   a   computer.      He    then   examined   the    simulated 
wave    records    for   surface   elevation , crest-to- trough  wave 
heights    and   zero-up-cross   wave   heights.      He    found   that    for 
the    zero-up-crossing   determination   of  wave   heights ,    the   Ray- 
leigh   distribution   is    a   good   approximation    irrespective   of 
the    spectral  width. 

Goda    (19  70)    suggests    that    the    deviation    of    the   wave 
height    distribution    from   the   Rayleigh      type    in   shallow  water 
should  be    considered   separately    from   the    change    in    the    func- 
tional  shape    of   the   wave    spectrum.       It   must   be    a  highly  non- 
linear phenomenon,    which   requires    full    analysis   with    the 
theory   of  non- linear  wave    spectra. 

Bretschneider    (1959)    investigated   the   PDFs    of  wave 
heights,    H,    wave    lengths,    L,    and  wave    periods,    T.       He    con- 
cluded  from  his    empirical    study    that    in    deep   water  wave 
systems,    the   PDFs    of  wave   heights    and    lengths    can   be    approxi- 
mated by    a   Rayleigh    distribution.       Using    the    relationship, 
wave    length   proportional    to   period   squared    for    deep    water 
linear   waves,    he    concluded   that    the    PDF    of  periods    squared 
is    a   Rayleigh    PDF   given   by 


p(T)    =    2.7 


exp 


■0.675 


4 


(3) 


where   f   is  the  mean  period. 
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Using   simulated  wave    data,    Goda    (1970)    found   that    the 
crest-to-crest   periods    squared  were    generally  not  well    ap- 
proximated by   the   Rayleigh   PDF.      The   mean    and   standard   devia- 
tions   of    the   periods   by    this    method  were    found    to   be    a 
function   of   the    spectral  width.      He    found   that    the  periods 
determined  by    the    zero-up-cross    method   are    only    in    fair 
agreement  with   Bretschneider    (1959).       If   the   waves    are   of 
the    type   which   have    sharply   peaked    spectra,    double    peaked 
spectra,    or   a   flat    spectrum,    Bretschneider ' s    (1959)    distri- 
bution was    totally   in    disagreement   with    that   of    the    simu- 
lated waves . 

Contrary    to   Bretschneider    (1959)    and   Goda    (1970), 
Longuet-Higgins     (1975)    noted   that    in    general    the    deep   water 
wave   periods   were   bell    shaped    (but  not   Gaussian)  .      He    defines 
the   period  PDF  when    calculated  by    the    zero-up-cross   method 
as  , 

1  3 

p(n)    =  y~  9    >      n  =    (t-<t>)/v<t>  (4) 

2(l+n    )    Z 

where      <t>      is    the   mean   period   and      \>      is    proportional    to 
the    spectral  width  where 

v   =    (u2/mo)%(<t>/2tt)  (5) 

and   u    and   |i«   are  the  zero  and  second  moments  of  the 
energy  spectrum. 

Koele  and  de  Bruyn  (1964),  Goda  (1967),  and  Siefert 
(1970)  have  observed  that  the  distribution  of  wave  heights 
in  shallow  water  does  not  correlate  well  with  the  Rayleigh 
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distribution.   Kuo  and  Kuo  (19  74)  have  suggested  that  this 
is  due  to:  1)  the  non-linear  effects  of  wave  interactions 
yielding  more  larger  waves,  2)  the  effect  of  bottom  friction 
yielding  a  reduction  in  the  low  frequency  components,  and 
3)  the  effects  of  wave  breaking,  which  would  truncate  the 
distribution  and  transfer  some  of  the  kinetic  energy  to  the 
high  frequency  components.   Therefore  the  PDFs  of  wave 
heights  with  a  certain  intensity,  may  be  considered  to  ap- 
proximate a  conditional  Rayleigh  distribution  truncated  by 
the  local  breaking  height. 

A  zero-up- cross  analysis  was  conducted  on  waves  measured 
on  a  reef  by  Black  (1978),  who  then  compared  the  observed 
distribution  with  the  Rayleigh,  truncated  Rayleigh,  and 
the  Weibull  distributions.   He  observed  that  the  truncated 
Rayleigh  distribution  fit  well  for  those  waves  in  the  breaker 
zone  when  the  heights  were  depth  limited.   The  truncated 
Rayleigh  distribution  did  not  fit  well  for  those  waves 
measured  either  inshore  or  offshore  of  the  breaker  zone. 
He  has  found  that  a  Weibull  distribution  more  closely  approx- 
imated these  measured  distributions. 

It  would  be  expected  that  the  Weibull  distribution  would 
yield  the  better  fit.   The  Rayleigh  distribution  is  a  func- 
tion of  the  variance  of  the  data,  whereas  the  Weibull  dis- 
tribution is  a  function  of  higher  moments  about  the  mean. 
Neither  theoretical  distribution  adequately  describes  ob- 
served distributions  which  are  multi-modal. 

Black's  (1973)  observations  of  wave  periods  showed  a 
distribution  which  was  skewed  to  the  right.   This  does  not 
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fit  well   with   the   symmetrical   Longuet-Higgins    (1975)    distri- 
bution.     It   was    found    that    the   period   distribution   more 
closely    followed   the   Rayleigh    distribution    than    the   period 
square    as    proposed  by    Bretschneider    (1959).       It   was    found 
that    the   Weibull    distribution   most    closely    fit    the   empiri- 
cal,PDF,    but    there    is    considerable   variation    in    the   peaked- 
ness    parameter.      His    attempts    in   relating    the   period  statis- 
tics   to    the   position   on    the    reef  were   unsuccessful.      He 
feels    that    this    is    due    to    the  non-linearities    of   the  waves 
in   shallow  water. 

In    this    thesis,    the   empirical    data   collected    for  breaking 
waves    were    compared  with    the    Rayleigh    distribution. 

D.       JOINT    PROBABILITY   DENSITY   FUNCTIONS 

A   study   of    the    joint    PDF    for  wave    heights,    n,    and    lengths, 
X,    which    can   be   viewed   as    heights    and  periods)    was    accom- 
plished by    Bretschneider    (1959).       The    general    form  of    the 
joint   PDF    for   two    related  variables    (r\ ,  X)    is    the   product   of 
the   marginal    PDF    (p(n))    and    the    conditional    PDF    (p(,\)). 

n 
p(n,A)  =  p(n)-p   (X)  (6) 


then 


00  00 

/  p(n,A)dA  =  P(n)     /  ; 


P(n)    -    /  p(n,A)dA   -  P(n)    /  pp(A)dA  (7) 

o  o 

Assuming  that  the  marginal  PDFs  are  Rayleigh,  Bret- 
schneider (1959)  examined  the  extreme  cases  of  0  and  +1  cor- 
relation.  For  the  case  of  zero  correlation,  the  marginal 
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Rayleigh  PDFs  are  given  by 

p(n)  =  I   n  exp(-nn2/4)  ;   n  =  H/S  (8) 


p(A)  =  y  X  exp(-TTA2/4)  ;   A  =  L/C  (9) 


For  zero  correlation,  i.e.  total  independence  (Figure  3), 
the  joint  PDF  is , 

p(n,A)  =  P(n)-p(A)  (10) 

so    that    the    joint    PDF    of   heights    and    lengths    is    given    as, 
p(n,A)    =   4   n  exp(-Trn   /4)  •  Aexp(-irA   /4)  (11) 

In    terms    of  heights    and   periods,    equation    15    is, 

p(n,t)    =    1.35  tt  exp(-Trn2/4)  •    t3   exp(-0.675    t4)  (12) 

where      t    =    T/f. 

For  the  case  of  total  dependence,  the  correlation  coef- 
ficient is  equal  to  one.   For  a  correlation  coefficient  of 
one,  all  data  points  on  a  plot  of  joint  Rayleigh  PDFs  fall 
on  a  45  degree  straight  line  passing  through  the  origin, 
n=A=0.   Assuming  that  a  joint  Rayleigh  distribution  is  appli- 
cable, the  following  relation  applies: 

p(n,A)  =  j   n  exp(-TTn  /4)  •  j   A  exp  (-ttA  /4)       (13) 

which  is  plotted  in  Figure  3. 
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Figure  3.   Statistically  Independent  Joint  Rayleigh  PDF. 
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Collins    (1970)    studied   the   PDFs    of  breaking  waves    de- 
rived  from   deep   water  waves   having   a  joint   Rayleigh    distri- 
bution  of  heights,    lengths    and   angles.      He    assumed   that 
H   =    yd    ,    where      II      is    the   breaker  height,       d      is    the    breaker 
depth,    and      y      is    the   breaking    criterion   assumed.      The    limit- 
ing  factor   of  waves    breaking   in    the    surf   zone   yields    a   dis- 
tribution  better   approximated  by   a   truncated  Rayleigh    dis- 
tribution.     As    expected,    he    found  that    the    longest    deep   water 
waves   were    associated  with    the    largest   breaking  waves. 

Collins     (1970)    concluded    from   the    study    of    the    joint 
PDFs    of  wave   heights    and    lengths,    that    the    effect    of    deep 
water    angle    of    approach    is    relatively   minor    for   breaking 
waves.       The   breaking  wave    steepness    was    found    to    be    very 
sensitive    to    the    choice    of   breaking    criteria. 

A  theoretical  expression  for  the  joint  distribution  of 
wave  periods  and  amplitudes  was  derived  by  Longuet-Higgins 
(19  75)    for   a  narrow   frequency   spectrum 

2 

p(C,n)    =  -£— c  exp(-C2d+n2)/2)  (14) 

(2tt)'2 

where      E,      and      n      are    the    normalized   amplitudes    and  periods 
as    defined  by, 

C    -    a/uo%  (15) 

and 

n  =    (t-<t>)/v<t>  (15) 

where      v      is    defined   in   equation    5.       The    wave    amplitude 
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is      a      ,       y        is    the    zero        moment   of   the   energy   spectrum, 
and      <t>      is    the   mean  wave   period.      Longuet-Higgins    (1975) 
compared   this    theoretical    distribution   of   ocean  waves    to 
that    of   the    deep   water   data   obtained  by   Brets chneider    (1959) 
and   found   good   agreement. 

Chakrabarti    and   Cooley    (19  77)    compared    the    theoretical 
joint   distribution   of  Longuet-Higgins    (1975)    with  wave    data 
recorded   during   a    1961   storm   in    the   North  Atlantic.      The 
spectrum  of   these   waves    did  not    fall    in    the    category   of 
narrow-band.      Through    the    comparison,    it   was    found   that 
there   was    a   considerable    agreement    at    the   higher  wave   heights 
but    a    definite    trend   away    from   the    theoretical    curve    at    the 
lower  wave   heights. 
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III.   EXPERIMENT 

A.   EXPERIMENTAL  SITES 

Three  different  California  experimental  sites  were  used 
to  obtain  data  on  the  three  types  of  breaking  waves --plung- 
ing, spilling,  and  collapsing.   During  the  period  4  March  to 
10  March  1975,  plunging  and  spilling  waves  were  measured  by 
Galvin  (19  75)  on  Del  Monte  Beach,  Monterey,  California. 
This  locale  was  chosen  because  the  waves  are  highly  refracted 
and  directionally  filtered,  yielding  narrowbanded  swell  off- 
shore which  impinges  almost  perpendicularly  on  the  beach. 

On  29  May  19  75  Galvin  (19  75)  conducted  a  second  experi- 
ment on  the  Carmel  River  Beach,  which  is  located  approximate- 
ly five  miles  south  of  Monterey,  California.   This  area  is 
in  an  embayment,  where  the  waves  are  again  highly  refracted 
and  directionally  filtered.   The  beach  at  this  site  is  very 
steep  resulting  in  plunging  and  often  collapsing  and  surging 
breakers.   Due  to  the  steepness  of  the  beach  and  the  rapidity 
of  breaking  of  the  shoaling  waves,  reflected  waves  are  found 
to  be  present. 

Spilling  breakers  on  a  gently  sloping  beach  were  measured 
in  March  1977  at  Torrey  Pines  Beach,  a  site  just  north  of 
La  Jolla,  California  (Olsen,  1977). 

On  3  and  9  March  19  78  data  were  collected  by  this  author 
on  Del  Monte  Beach,  Monterey,  California.   On  the  first  day 
of  the  experiment,  very  narrowbanded  swell  was  observed 


24 


breaking   on   the   shore,    mainly   of    the   plunging    type.      A   storm 
arrived   late   on   8  March   and   continued  through    9   March.      Winds 
of  approximately   15   knots    yielded  a   very    confused   sea. 
Breakers    were    of  both    the    spilling    and   plunging    type    and 
were    observed   to   break    at   various    small    angles    relative    to 
the  beach.       Sea   surface   elevations    and   velocities   were    col- 
lected  for   approximately    three  hours    on   each   of   these    days 
on   both    the    flowing   and  ebbing   sides    of  high    tide. 

B.       INSTRUMENTS 

Capacitance    type    gauges   were   used   to   measure    the    change 
in    the   sea   surface    elevation    as    explained  by   Thornton   et    al. 
(1976).       The   wave    gauges    were    constructed   of    3/8    inch    poly- 
propylene   covered   stainless    steel    rod.       Calibration  was 
accomplished   in    the    laboratory    statically   by    lowering    the 
wave    gauge    into    a    deep   water    tank.      Depth    of    immersion   was 
then   plotted   against    volts. 

The   velocity   of   the   water  particles   was    measured  with 
Marsh-McBirney  Model    721    and    722   Electromagnetic   Current 
meters.       These    flowmeters    work   on    Faraday's    principle    of 
electromagnetic    induction    and    are    capable    of    recording    the 
particle   velocities    in   orthogonal    directions    in    a  plane. 
The   meters    were    dynamically    calibrated   in    the    laboratory. 
Calibration  was    accomplished  by   oscillating    the    current 
meters    in    a  water    tank  with    a  prescribed  motion    simulating 
the    flow   conditions    which   would  be    expected    to    be    found    in 
the    surf   zone.       The    flowmeter    voltage    output   was    found    to 
be    linear    in    the    velocity    ranges    found   in    the    breakers. 
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All    instruments   were   placed   on    towers    in    the    surf   zone 
at    low   tide    at    the   expected  point   of  breaking   of   the   waves 
for   the    subsequent   high    tide.      Each   tower   held   a    capacitance 
wave   gauge    and   an   electromagnetic   flowmeter.      The    flowmeter 
was    placed  below   the  wave    gauge,    approximately    0.5   m   from 
the   bottom,    to    insure    constant    submergence.      Data  were    col- 
lected within   an   hour   or   two   on   both    the   ebb    and   flow   sides 
of  high    tide.      Analog   recordings    were    made    on   both    a   strip 
chart   and  magnetic   tape. 

C.       DATA  ANALYSIS 

All   signals   were    digitized   from   their   analog    form  by 
means    of   a  Vidar   Corporation    32-channel    digital    data   acqui- 
sition   system.       The    digitized    data  were    then   processed   on 
an    IBM   360/67    computer.      Record    lengths    of   approximately 
30   minutes    from  each    data   set   were    analyzed.       Specific   de- 
tails   of    the    analysis    can   be    found   in   Appendix  A. 

The    data  were    edited  by    excluding    any    value    greater 
than    six   standard    deviations    from    the    mean    and   replacing 
that    value  with   the   previous    value.      The   values   were    then 
linearly    detrended    to    exclude    the    effects    of    the    rising    or 
falling   tides.      The   mean    and  variance   of   sea   surface    eleva- 
tions  were    calculated    for   each    record   as    well    as    the    mean 
and  variance    of    the   heights,    periods    and   maximum  velocities 
all   of  which   are    defined   on   page  27. 

It   was    desired    to    only    examine    the    sea-swell    band   of 
frequencies    lying   between    0.03-1.0    Hz    (30    sec-1    sec).    Hence, 
the    data  were    passed   through    both    a   high    and    low   pass    digital 
filter. 
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The  high  pass    filter   had   a   cut-off   frequency   of   0.03   Hz 
(30    sec).      The   high  pass    filter   used   a   Fast   Fourier  Trans- 
form   (FFT)    algorithm   to   obtain   the   Fourier   spectrum  of   the 
entire    30   minute    record.      The   Fourier    coefficients    correspond- 
ing  to   0-0.03   Hz   were    then   used   to   synthesize    a    low   frequency 
time    series   which  was    subtracted   from   the  wave    record. 

The  low  pass  filter  was  a  25  weight  inverse  transform 
filter  (Davidson,  1970),  as  shown  in  Figure  4.  A  cut-off 
frequency  of  1.2  Hz  was  used  with  a  terminal  frequency  of 
1.8    Hz. 

Wave   heights   were    determined   from   the    surface    elevation 
record  by   means    of    the    zero-up-cross    method.       The    zero-up- 
cross    technique    defines    a  wave   period   as    the    interval   between 
adjacent   upcrosses.      The   maximum   and  minimum  of   the    surface 
elevation   within    the    interval    defines    the    crest    and    trough 
of   a  wave.      Maximum  onshore   particle   velocities   were    calcu- 
lated   from   the    flowmeter   records.       Maximum   onshore    velocity 
was    defined   as    the   maximum  velocity   amplitude    occurring  be- 
tween  successive    zero   upcrossings    as    defined  by    the    sea 
surface    elevation    record.       Joint    PDFs    of  wave    height    and 
period,    wave   height    and  maximum  velocity,    and  maximum  veloc- 
ity   and   period  were    calculated    from   the    analyses. 

Heights    and  maximum  velocities    were   normalized  by    divid- 
ing  by    their    respective    standard    deviations.       The    periods 
were   normalized  by    dividing   by    their   mean    period.       The   nota- 
tion   us e d    is 

H  -   H'/cr 

n 

V  =   V»/av  (17) 
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Figure  4.   Low  Pass  Digital  Filter, 


where   heights    (H1  )  ,    maximum  velocities    (V),    and  periods 
(T1)    are    the   values    as    determined  by    the    zero-up-cross 
method.      f      is    the   mean   period,    and      a        and      a        are    the 
standard   deviations    of   the    surface    elevations    and  particle 
velocities    respectively.      Histograms   were    calculated   for 
the  normalized  heights,    periods,    and  maximum  velocities. 
The   histograms   were    transformed   into    PDFs    by    dividing    the 
number   of  value   over   an    interval  by    the    total   number   of 
values    in   the    array   and  by    dividing  by    the    interval  width. 
The  mathematical   expression    for   the   exact   PDF   is 

rs 


p(x)    =    lim   x^Q    J_    {li 


lm 


Ax  N-*»  |_N  _ 


x 


}  (18) 


where   N    is  the  number  of  times  the  variable  occupies  an 
interval,   N   is  the  total  number  in  the  record  observed, 
Ax   is  the  interval  width  and  p (x)   is  the  PDF.   A  restric- 
tion on  the  probability  density  is  that 


/  p(x)dx  =  1  . 
o 

Joint  probability  density  functions  of  periods  versus 
heights,  velocities  versus  heights  and  periods  versus 
velocities  were  calculated  in  a  similar  fashion.   A  15  by  15 
grid  was  used  to  determine  the  frequency  of  values  which 
simultaneously  fill  each  joint  interval.   Mathematically 
the  probability  that  the  arrays   x   and   y   simultaneously 
assume  values  within  the  range  (x,  x+Ax) ,  (y,  y+Ay)   may  be 
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obtained  by  taking  the  ratio  of  the  total  number  of  times  of 
each  occurrence  within  a  bin,  N   ,  and  the  total  number  of 
observations ,  N  .   The  exact  probability  density  is  approached 
as  the  number  of  observations  approaches  infinity  the  interval 
widths  approach  zero.   The  joint  PDF  can  be  defined  as: 


p(x,y)  -  lim 


1   rlim 

x  -*■  °°   AxAy    N 

y  -»■  co 
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2SX 
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(19) 


with  the  restriction 

-f-co         -|-co 


// 


O  O 


p (x,y)    dx    dy   =    1 


The   energy    density    spectrum    (power   spectrum)    was    obtained 
by   means    of   a   Fast    Fourier   Transform    (FFT)    algorithm.       The 
energy    density   spectrum   tells    how   the   variance    is    distributed 
with   respect    to    frequency. 

The    spectral  width   parameter  was    calculated   using    the 
equation: 


e    = 


Jm0m4" 
f  m0m4 


a0m4"m2 


(20) 


where      m~  ,    m« ,    and  m,       are    the    spectral    moments     (Cartwright 
and  Longuet-Higgins ,     1956).       For    e=0    the   waves    have    only    one 
frequency    and    the    energy    spectrum  would  be    represented   by    a 
spike.      As    £-*l    the    sea    surface    can    be    described    as    a   broad- 
band process    for  which    the    distribution   would   be    Gaussian 
and   the    spectrum  would   approach    a   horizontal    line.       Swell 
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can   therefore   be    described   as    a  narrow-band  process    and   seas 
as   broad-band.      Due    to    the   appearance    of   secondary  waves ,    a 
narrow-banded   deepwater   spectrum  of  waves   broadens    in    the 
surf   region. 

The    correlation    coefficients    of   the   joint   values   were 
calculated   to   obtain   a  measure    of    the   amount    of    linear   rela- 
tionship   of  random  pairs    of  values    of  periods    and  heights  , 
velocities    and  heights,    and  periods    and   velocities, 

r(x.7)    =   NTT    D  a        a  ^ 

1-1  x        x 

where      x      and     y      are    the    sample   means,      x.       and     y.       are 
the    sample    standard   deviations. 
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IV.      DISCUSSION 

Miller    (1976)    concluded   that   breaker   shapes    form   a    con- 
tinuous   gradation   of   shapes    from  spilling    to   plunging.       Each 
breaker   type   undergoes    the   same   processes,    but    differs    only 
in   the  magnitude    of   the   vortices    generated.      The   probability 
density    functions    appear   to    show   a   similar   gradation.         Two 
simplified   cases,    for   spilling   and  plunging,    were    chosen    to 
represent    extreme   examples    and   are    used   for   all    discussions. 
The    4  March    1975    data   most    simply    represent    the   plunging 
case    and   the    16   March    1977    data   represent    the    spilling    case. 
The    results    for    all    data    are    given    in   Appendix  C. 

In   both    examples    for    spilling    and   plunging   breakers,    the 
offshore   waves    appeared   as   narrow   band   swell.      A   spectral 
analysis    of   the    surface   elevation    records    in    the    surf   zone 
yielded   spectral  width   parameters    ranging    from  0.94-0.97 
(see   Appendix  A),    suggesting    a   broad  band   spectrum   in    the 
surf   zone.      The   broad   spectral  width    can   be    attributed   to 
the   non- linear  wave-wave    interactions    yielding   harmonics    in 
the   breaking    zone. 

In   order    to    determine    the    origin   of    the    large    spectral 
width  parameter,       e      was    calculated   in    the    ranges    of   0-0.25 
Hz,    0-0.50    Hz,    0-1.0    Hz    and   0-nyquist    frequency.       The    analy- 
sis   revealed   that    large   values    of   the    spectral  width   para- 
meter   occurred   in    all    cases.       In    the    case    of   0-0.25    Hz   which 
contains    the    first    three    harmonics,       e      was    greater    than    0.7, 
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which  was    the    case  with    the    smallest      e    .       It    is    concluded 
that    the    change    from  a  narrow   deep   water   spectrum   to    a 
broad-band   shoaling   spectrum  is   primarily    due    to    generation 
of   secondary  waves    at   harmonic    frequencies. 

A.       PERIOD    PROBABILITY  DENSITY   FUNCTIONS 

Bretschneider    (1959)    made    the    assumption    that    the   wave 
lengths   were   Pv.ayleigh    distributed;    using    this    assumption 
and   the    linear  wave    theory   relationship 

L  =   £  T2  (22) 

■ 

the  wave  periods  squared  would  be  Rayleigh  distributed  in 
deep  water.   Bretschneider' s  analysis  can  be  extended  to 
the  shallow  water  case  using  the  shallow  water  approxima- 
tion for  wave  length 

L  =  >^h  T  (23) 

Assi.Tmi.ng  the  wave  lengths  are  Rayleigh  distributed,  the 
periods  in  shallow  water  would  be  expected  to  be  Rayleigh 
distributed. 

The  period  PDFs  for  plunging  and  spilling  examples  are 
shown  in  Figures  5  and  6  respectively.   The  Rayleigh  PDF  has 
been  superimposed  on  the  period  PDFs  for  comparison.   The 
formulas  used  for  the  Rayleigh   PDF  for  various  parameter 
transformations  are  given  in  Appendix  A.   In  both  the  spill- 
ing and  plunging  wave  examples,  it  is  noted  that  the  period 
PDFs  only  loosely  follow  the  Rayleigh  PDF.   The  Rayleigh  PDF 
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overpredicts  the  wave  periods  greater  than  the-  Rayleigh  mode 
value,  and  underpredicts  the  observed  periods  less  than  the 
peak  Rayleigh  value.   The  truncation  of  the  long  period 
waves  is  probably  due  to  the  larger  period  waves  breaking 
farther  offshore  and  then  reforming  to  be  recorded  by  the 
stationary  wave  gauge  as  higher  frequency  waves. 

The  spilling  example  for  16  March  19  77,  shown  in  Figure  5, 
has  the  mode  lying  close  to  the  maximum  of  the  superimposed 
Rayleigh  PDF.   The  peak  of  the  breaking  wave  distribution 
is  located  at  the  first  harmonic  of  the  offshore  wave  period, 
which  is  attributed  to  the  secondary  waves.   A  peak  corre- 
sponding to  the  peak  offshore  wave  period  is  absent. 

The  plunging  waves,  as  represented  by  the  4  March  19  75 
record  in  Figure  6,  are  multimodal.   Since  the  Gaussian,  Gram- 
Char  lier,  Rayleigh  and  Weibull  PDFs  are  unimodal,  none  of  the 
theoretical  PDFs  will  adequately  describe  multi-modal  PDFs 
for  this  particular  plunging  surf  condition. 

Two  strong  peaks  are  found  to  occur  in  the  plunging 
breaker  period  PDF.   The  longer  period  is  found  to  correlate 
well  with  peak  period  of  the  wave  spectrum.   The  shorter  wave 
period  appears  at  approximately  the  first  harmonic.   Since 
the  center  period  of  the  bin  is  used  to  represent  the  inter- 
val, the  magnitude  of  the  peaks  and  the  proximity  to  har- 
monics as  represented  by  the  labels  is  a  function  of  the 
interval  width  and  the  interval  cutoffs.   A  30  bin  PDF  was 
calculated  for  better  resolution  to  determine  if  the  peaks 
were  harmonics  and  Figure  7  shows  that  they  are. 
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The  waves   periods    as    determined  by   a   zero-up-cross    of 
the   surface   elevation    record   for   spilling   conditions    yields 
a  poor   indication   of    the   peak   of   the   energy   spectrum  which 
corresponds    to    the   mean   offshore   wave   period.       In    the    spill- 
ing breaker   case,    the   mean  period   corresponded   to    the    first 
harmonic   of  peak    in    the    energy    spectrum.      For    the   plunging 
breakers,    the   mean   period    corresponded   to    a  valley   in   the 
period  PDF  between    the   primary   and   first   harmonic  periods. 
The    appearance    of   strong   peaks    at   one   half   the   offshore   wave 
period   is    indicative   of   the    importance    of  non-linearities 
in   the    form  of   secondary  waves    in    the    surf   zone. 

B.       HEIGHT   PROBABILITY   DENSITY   FUNCTIONS 

Wave   height    PDFs    for   the   plunging    and   spilling  wave    ex- 
amples   are    shown   plotted   against    the    theoretical   Rayleigh 
PDF    in    Figures    5    and   6.       The   height   PDF    for   spilling  waves 
is   not    as    easily    described   as    the    period   PDF.       The    height 
PDF    is    loosely    described  by    the    Rayleigh   PDF.       The    major 
height    peak    is    centered   at    2.2.       The    Rayleigh   PDF    overesti- 
mates   the    density    of    the    small  wave    observations    and    tends 
to    underestimate    the    density   of    the    large   waves.       The   heights 
are    truncated   at    approximately    six   standard    deviations.    Since 
the   wave    gauge    is    stationary    and    therefore    the    depth   of    the 
water   is    essentially    fixed,     this    truncation    occurs    due    to 
the   breaking   of    the    large   waves    seaward   of    the   wave    gauge. 
The    calculated  mean    height    value    on    16    March    19  77    lies    at 
2.7,    which    actually    corresponds    to    a    dip    in    the    PDF    relative 
to    the   peak    intervals    on    either    side. 
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The  mean  height  observed  in  the  4  March  19  75  plunging 
breaker  PDF  was  calculated  to  be  at  2.4.   The  observed  PDF 
is  found  to  be  strongly  bimodal  with  peaks  at  1.8  and  3.4. 
The  strong  bimodality  is  attributed  to  the  very  narrow  band 
incident  wave  system  and  the  generation  of  secondary  waves 
at  the  first  harmonic. 

C.   JOINT  PROBABILITY  DENSITY  FUNCTIONS  (T  vs  H) 

The  joint  PDF  of  periods  and  heights  for  the  plunging 
and  spilling  examples  are  contoured  in  Figures  5  and  6;  the 
distributions  for  all  cases  are  found  in  Appendix  C. 

It  is  noted  that  the  longest  wave  periods  are  associated 
with  the  greatest  wave  heights  suggesting  strong  correlation 
between  periods  and  heights.   The  correlation  coefficient 
was  calculated  for  all  cases  and  varied  between  0.60  and  0.80, 
Bretschneider  (1959)  examined  the  two  special  cases  of  joint 
Rayleigh  distributions  with  correlations  of  0  and  +1.   Per- 
fect correlation  between  two  variables  would  result  in 
correlation  coefficient  of  1.0  and  a  straight  line  in  the 
joint  PDF. 

The  joint  Rayleigh  PDF  for  the  case  of  a  correlation  co- 
efficient of  zero  is  plotted  in  Figure  3.   Zero  correlation, 
or  statistical  independence,  results  in  the  joint  PDF  being 
spread  out  as  compared  with  the  straight  line  for  perfect 
correlation.   The  amount  of  spread  is  a  measure  of  the  reduc- 
tion in  correlation. 

Figure  5  is  an  example  of  p(H,T)  for  plunging  breakers 
and  shows  two  regions  of  high  density.   The  bimodality  of 
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p(H)  and  p (T)  are  reflected  in  the  p(H,T).   The  maximum  peak 
period,  which  corresponds  to  the  peak  in  the  wave  spectrum, 
is  highly  correlated  with  the  primary  peak  of  wave  heights. 
The  second  high  density  region  corresponds  to  the  first 
harmonic  of  the  primary  wave  period.   The  joint  PDF  suggests 
strong  wave-wave  interactions  in  a  plunging  wave  system  re- 
sulting in  energy  transfer  to  harmonics. 

Two  distinct  contour  groups  are  no  longer  apparent  in 
Figure  6  for  the  spilling  case,  but  rather  one  long  contour 
is  found.   Since  the  height  PDF  and  the  period  PDF  are  only 
weakly  bimodal,  it  follows  that  a  single  contour  group 
should  be  found.   The  area  of  maximum  joint  density  is  found 
to  be  associated  with  the  first  harmonic  of  the  offshore 
wave  period,  and  the  mode  of  the  observed  breaker  heights. 

Greater  density  was  found  at  short  periods  and  heights 
for  all  records  than  predicted  by  the  Rayleigh  PDF.   A  high 
frequency  wave  riding  on  a  long  period  wave  causes  an  addi- 
tional zero  crossing  as  the  longer  wave  passes  through  zero. 
Most  previous  comparisons  of  heights  or  periods  of  waves  have 
been  done  using  pressure  sensors  which  greatly  filter  the 
high  frequencies  due  to  hydrodynamic  attenuation.   Capaci- 
tance wave  staffs  and  electromagnetic  flowmeters  have  good 
response  times  and  thus  measure  high  frequencies.   Therefore, 
even  though  the  data  were  filtered  at  1  Hz,  all  records  showed 
peaks  at  the  shortest  periods  and  heights. 
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D.  VELOCITY  PROBABILITY  DENSITY  FUNCTIONS 

The  velocity  PDFs  were  calculated  as  the  maximum  onshore 
velocities  within  the  period  segments  determined  by  the  zero- 
up-crosses  of  the  surface  elevation  record  as  shown  in  Fig- 
ures 8  and  9.   Using  this  technique  can  result  in  both  posi- 
tive and  negative  maximas  for  onshore  particle  velocities 
(i.e.  ,  d2v/dt  <  0). 

The  velocity  PDFs  for  spilling  breakers  and  plunging 
breakers  are  shown  in  Figures  8  and  9.   The  theoretical 
Rayleigh  PDF  is  superimposed  on  the  same  graphs  for  compari- 
son.  It  is  apparent  that  neither  the  spilling  or  plunging 
case  is  adequately  described  by  the  Rayleigh  PDF.   The  veloc- 
ity PDF  for  spilling  breakers  appears  to  be  symmetrical,  with 
the  mode  value  falling  at  the  bin  with  center  velocity  of 
1.5.   There  also  appears  to  be  a  small  peak  value  at  0.3. 
When  compared  with  the  plunging  case,  it  is  found  that  the 
peaks  occur  at  the  same  normalized  velocities.   The  plunging 
wave  PDF  is  bimodal  and  has  relatively  few  occurrences  of 
negative  onshore  velocities. 

E.  JOINT  PROBABILITY  DENSITY  FUNCTIONS  (V  vs  H) 

The  joint  probability  distributions  of  velocities  and 
heights  are  shown  in  Figures  8  and  9.   The  results  show  high 
correlation  between  velocities  and  heights.   The  high  posi- 
tive correlation  suggests  that  the  larger  maximum  onshore 
velocities  are  correlated  with  the  larger  wave  heights. 

Figure  3  shows  the  plunging  breaker  case.   Two  peak  areas 
are  found  to  occur  in  the  joint  PDF  which  reflects  the 
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bimodality  of  the  individual  PDFs .   The  joint  PDF  for  the 
plunging  waves  suggests  a  highly  correlated  case  (0.68)  in 
which  the  largest  onshore  velocities  occur  due  to  waves  of 
greatest  height. 

The  spilling  case  as  shown  in  Figure  9  also  shows  high 
correlation  (0.70).   For  this  example,  there  is  only  one 
major  peak  which  corresponds  to  the  individual  PDFs  which 
are  essentially  unimodal.   The  peak  area  of  the  joint  PDF 
has  been  shifted  toward  higher  waves  and  greater  velocities 
than  the  plunging  case. 

F.   JOINT  PROBABILITY  DENSITY  FUNCTIONS  (T  vs  V) 

Figures  10  and  11  show  joint  PDFs  of  periods  and  veloc- 
ities in  which  the  periods  were  determined  from  the  zero-up- 
cross  es  of  the  wave  gauge  record  and  the  coincident  maximum 
onshore  velocities  are  determined  from  the  flowmeter  record. 
The  results  are  similar  to  the  joint  PDF  of  period  and 
heights.   Spilling  waves,  as  contoured  in  Figure  11,  show  one 
main  contour  group.   The  maximum  spilling  probability  is 
found  at  the  intersection  of  period  0.84  and  velocity  1.5. 
The  plunging  waves  as  found  in  Figure  10  reveal  two  contour 
peaks.   One  large  contour  group  is  associated  with  the  point 
representing  the  offshore  period  1.56  and  a  velocity  of  1.5. 
The  second  large  contour  is  a  function  of  the  first  harmonic 
wave  period  of  0.84  and  a  velocity  of  0.3. 

The  negative  values  of  the  maximum  velocity  indicate 
that,  during  a  wave  period  defined  by  the  zero-up-crosses 
of  the  surface  elevation,  the  maximum  velocity  can  be  negative, 
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The  joint  PDF  of  wave  height  and  velocity  shows  that  the 
negative  velocities  are  associated  with  the  low  wave  heights 
The  joint  PDF  of  velocity  and  period  shows  that  small  nega- 
tive velocities  are  correlated  with  the  short  periods . 
Therefore,  it  is  concluded  the  negative  velocity  maxima 
are  associated  with  short  periods  and  low  wave  heights.  The 
secondary  waves  occur  after  the  passage  of  the  primary  wave 
crest  and  have  a  peak  near  the  trough  of  the  primary  wave; 
the  velocities  associated  with  the  primary  wave  would  be 
negative  maximum  in  this  region  so  that  the  measured  veloc- 
ity, which  is  a  sum  of  the  primary  and  secondary  waves,  can 
be  small  or  negative. 
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V.   CONCLUSIONS 

Waves  in  the  surf  zone  are  a  highly  non-linear  process 
which  is  evidenced  by  the  appearance  of  secondary  waves. 
The  secondary  waves  appeared  in  the  PDFs  at  the  first  har- 
monic of  the  primary  period. 

The  joint  probability  density  functions  of  the  breaking 
wave  properties  show  high  correlations.   The  correlation 
factor  ranged  between  0.60-0.80  for  periods  and  heights 
and  suggests  that  the  greater  the  wave  period,  the  higher 
the  breaker  height  will  be.   High  correlation  between  the 
periods  and  velocities  and  the  heights  and  velocities  sug- 
gests that  the  maximum  onshore  velocities  are  a  function 
of  wave  periods  and  breaker  heights. 

Breaking  waves  may  be  viewed  as  a  gradation  between 
two  extremes,  spilling,  and  plunging  breakers.   The  differ- 
ence between  the  spilling  and  plunging  breakers  becomes 
apparent  in  the  joint  PDFs.   The  spilling  breakers  yield  a 
joint  PDF  which  has  only  one  high  density  contour  region, 
whereas  the  plunging  breakers  show  two  high  density  contour 
groups.   The  high  density  contour  groups  appear  to  be  most 
associated  with  the  period  PDF. 

The  joint  PDF  (p(T,H))  of  spilling  breakers  has  a  maxi- 
mum density  at  a  wave  period  equal  to  one  half  the  peak 
period  of  the  wave  energy  density  spectrum,  i.e.  the  first 
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harmonic.   The  corresponding  velocities  (p(v,T))  or  heights 
(p(H,T))  were  likewise  found  to  be  unimodal  with  a  peak  den- 
sity which  is  driven  by  the  first  harmonic  wave  period. 

Plunging  breakers,  at  the  opposite  end  of  the  gradation, 
yield  a  bimodal  joint  PDF.   The  two  major  contour  groups 
appear  to  be  a  result  of  the  bimodal  density  distribution 
of  the  periods.   The  two  peak  densities  of  the  period  PDF 
are  representative  of  the  peak  of  the  wave  density  spectrum 
and  a  first  harmonic.   Due  to  the  high  correlation  of  periods 
and  heights ,  the  corresponding  PDF  of  wave  heights  appears 
to  be  bimodal.   The  joint  PDF  of  period  and  velocity  and  the 
joint  PDF  of  velocity  and  height  suggests  high  correlation 
between  these  variables.   The  maximum  onshore  velocities 
appear  to  be  correlated  with  both  wave  periods  and  wave 
heights . 

A  truncated  R.ayleigh  PDF  was  found  to  loosely  approxi- 
mate the  empirical  period  and  the  height  PDFs  of  spilling 
breakers.   This  truncation  seems  to  be  the  result  of  long 
period  waves  with  their  associated  larger  breaker  heights 
breaking  prior  to  the  arrival  at  the  wave  gauge.   The  mode 
of  the  spilling  wave  period  distributions  corresponds  to  the 
first  harmonic  of  the  peak  of  the  energy  density  spectrum. 
The  smoothed  appearance  of  the  spilling  wave  periods  PDF, 
relative  to  the  plunging  case,  appears  due  to  the  smearing 
which  occurs  with  the  appearance  of  harmonics  greater  than 
the  first.   The  velocities  PDFs  of  spilling  waves,  where  the 
velocity  maximum  is  determined  as  in  Section  D,  yields  a 
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symmetrical  distribution  with  some  negative  onshore  veloc- 
ities present.   Since  the  Rayleigh  distribution  is  undefined 
for  negative  values,  the  resulting  fit  of  observed  data  with 
the  theoretical  Rayleigh  distribution  is  poor. 

The  PDFs  of  the  plunging  breakers  are  poorly  described 
by  either  the  Rayleigh  PDF  or  other  existing  theoretical 
PDFs  such  as  the  Weibull  PDF.   Theoretical  PDFs  are  unimodal, 
whereas  the  empirical  PDFs  under  plunging  conditions  are 
strongly  bimodal.   The  high  correlation  found  in  the  joint 
PDFs  for  wave  heights  and  maximum  velocities  also  results 
in  a  bimodal  distribution  both  in  the  velocities  and  the 
heights  PDFs. 

The  mean  wave  period  as  obtained  by  the  zero-up-cross 
method  in  the  surf  zone  is  a  poor  indicator  of  the  mean 
offshore  wave  period.   Spectral  analysis  of  the  breaking 
waves  shows  an  energy  peak  corresponding  to  the  primary  off- 
shore wave  period,  but  the  period  PDFs  of  breaking  waves 
have  a  strong  peak  at  the  first  harmonic  of  the  wave  spectra 
peak.   Plunging  waves  had  high  density  peaks  in  the  period 
PDF  corresponding  both  to  the  peak  of  the  wave  spectrum  and 
to  its  first  harmonic.   The  result  is  that  the  mean  of  the 
breaking  waves  falls  between  the  primary  and  the  first  har- 
monic wave  period.   The  mean  period  of  breaking  waves  there- 
fore is  a  poor  estimation  of  the  mean  period  of  the  incident 
waves . 
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APPENDIX  A 
STATISTICAL  PARAMETERS  CALCULATED 

1.   Due  to  the  variation  of  sampling  frequencies  between  the 
different  experiments,  it  follows  that  some  of  the  other 
statistical  parameters  will  also  vary.   A  summary  of  these 
parameters  follows  where   f„   is  the  Nyquist  frequency,   Af 
is  the  resolution,   At   is  the  interval  between  samples  and 
T   is  the  length  of  sampled  record  in  minutes. 


DATE 

fN(Hz) 

Af    (Hz) 

At  (sample/sec) 

T(min) 

4  MARCH 

1975 

2.5 

6.1xl0~4 

0.200 

27.3 

5    MARCH 

1975 

2.5 

6.1xl0~4 

0.200 

27.3 

6   MARCH 

19  75 

2.5 

6.1xl0~4 

0.200 

27.3 

3  MARCH 

1975 

2.0 

4.8xl0~4 

0.256 

35.0 

29   MAY 

19  75 

2.0 

4.9xl0"4 

0.256 

35.0 

16   MARCH 

1977 

2.0 

4.9xl0"4 

0.250 

34.1 

17   MARCH 

1977 

2.0 

4.9xl0"4 

0.250 

34.1 

19   MARCH 

1977 

2.0 

4.9xl0~4 

0.250 

34.1 

3   MARCH 

1978 

2.0 

4.9xl0"4 

0.250 

34.1 

9   MARCH 

1978 

2.0 

4.9xl0"4 

0.250 

34.1 

Tables  of  calculated  statistical  parameters  follow  where 

H   =  sea  surface  elevation  (m) 

v  =  instantaneous  particle  velocities  (m/sec) 

VM  =  maximum  onshore  particle  velocities  (m/sec) 

LR  =  length  of  record  for  calculated  statistics  (min) 

e  =  spectral  width  parameter 

H   =  wave  height  (m) 

T   =  wave  period  (sec) 

r  =  correlation  coefficient 
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APPENDIX  B 
MATHEMATICAL  DERIVATION  OF  NORMALIZED  RAYLEIGH  DISTRIBUTIONS 

Wave  heights,  periods  and  velocities  have  been  compared 
with  the  Rayleigh  distribution 

A       ,    A2 


p(A)  =  -^  exp  (-  — j) 


4a         2a 

where      A      is    the  non-normalized  value.       Since    all   parameters 

are   normalized  in    this   paper,    the    following    conversion  was 

made : 

let      A'    =  A/a    , 

then  ^ 

p(A')    =   p(A)     (ig.)      A=A'a 


p(A')    =   p(A)    d(A\a)    =   p(A)a 


^A1 
for  wave   heights 

p(H)    =   -5-,    expC-A;^)2) 


4a    c 


2a 


p(^-)=zA-     exp(-%(2§-)2) 

n  n  n 


for  maximum  velocities 


for  periods 


p(V)  =  -X-2  exp(-%(2|-)2) 
4av  v 

V       V  V 


p(T)  =  -^  expC-i;^)2) 
4a 
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where  a  =  T/2 


then 


p(-)    =   -1?      exp(-%(^-)2) 
T  4a  ^a 


p(l)    =      I        exp(-%(^)2) 

m  qi  rn 
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APPENDIX  C 
PROBABILITY  DENSITY  FUNCTIONS 

The  following  tables  consist  of  both  the  normalized 
joint  probability  density  functions,  and  the  individual 
probability  density  functions  as  labeled.   The  PDF  tables 
coincide  with  Figures  5  thru  11,  respectively. 

The  top  table  is  the  joint  PDF.   The  top  line  and  the 
left  column  represent  the  center  period,  height,  or 
velocity  point  of  each  of  the  15  bins.   The  numbers  in 
the  other  columns  and  lines  represent  the  frequency  of 
occurrence  in  the  respective  bins. 

The  bottom  table  represents  the  individual  PDFs  as 
labeled.   The  top  line  of  the  individual  PDFs  is  the 
center  point  of  the  individual  bins.   The  second  line 
represents  the  frequency  of  occurrence  in  the  bin. 
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